Giant Thermoelectric Effect from Transmission Supernodes 
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We predict an enormous order-dependent quantum enhancement of thermoelectric effects in the 
vicinity of a higher-order 'supernode' in the transmission spectrum of a nanoscale junction. Single- 
molecule junctions based on 3,3'-biphenyl and polyphenyl ether (PPE) are investigated in detail. The 
nonequilibrium thermodynamic efficiency and power output of a thermoelectric heat engine based 
on a 1,3-benzene junction are calculated using many-body theory, and compared to the predictions 
of the figure-of-merit ZT. 
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Thermoelectric (TE) devices are highly desirable since 
they can directly convert between thermal and electrical 
energy. Electrical power can be supplied to such a device 
to either heat or cool adjoining reservoirs (Peltier effect) 
or alternatively, the flow of heat (e.g. from a factory 
or car exhaust) can be converted into usable electrical 
power (Seebeck effect). Often, the efficiency of a TE de- 
vice is characterized by the dimensionless figure-of-merit 
ZT=S 2 GT / k, constructed with the rationale that an ef- 
ficient TE device should simultaneously: maximize the 
electrical conductance G so that current can flow without 
much Joule heating, minimize the thermal conductance 
K in order to maintain a temperature gradient across the 
device, and maximize the Seebeck coefficient S to ensure 
that the coupling between the electronic and thermal cur- 
rents is as large as possible*^ Generally, however, ZT is 
difficult to maximize because these properties are highly 
correlated with one another,— ~— a fact that becomes more 
pronounced at the nanoscale where the number of degrees 
of freedom available is small. 

If a TE material were found exhibiting ZT>A it would 
constitute a commercially viable solution for many heat- 
ing and cooling problems at both the macro- and nano- 
scales, with no operational carbon footprint^ Currently, 
the best TE materials available in the laboratory ex- 
hibit ZTi=s3, whereas for commercially available TE de- 
vices ZT&1, owing to various packaging and fabrication 
challenges 

In a previous article, enhanced thermoelectric effects 
were found in the vicinity of a transmission node of a 
quantum tunneling device. Generically, the transmission 
probability vanishes quadratically as a function of energy 
at such a transmission nodeZ Here we present results for 
a class of two-terminal single-molecule junctions (SMJ) 
with higher-order 'supernodes' in their transmission spec- 
tra. In the vicinity of a 2n th order supernode: 



T(E) <x(E- /inode) 2 



(1) 



where (Anode is the energy of the node. We find that junc- 
tions possessing such supernodes exhibit a scalable order- 
dependent quantum-enhanced thermoelectric response. 

As an example, ZT of a supernode-possessing 
polyphenyl ether (PPE)-based SMJ is shown as a func- 



tion of repeated phenyl unit number n in Fig. {J}. As 
illustrated in the figure, ZT pea ^ scales super-linearly in 
n whereby ZT pea ^=A.l in a junction composed of just 
four phenyl groups (n=4). Although we focus on molec- 
ular junctions in this article, it should be stressed that 
our results are applicable to any device with transmission 
nodes arising from coherent electronic transport. 

As an engineering rule-of-thumb, ZT has been widely 
used to characterize the bulk thermoelectric response of 
materials At the nanoscale, however, it is unclear 
the extent to which ZT is applicable, since bulk scaling 
relations for transport may break down due to quantum 
effects.- Moreover, ZT is a linear response metric, and 
cannot a priori predict nonequilibrium thermoelectric re- 
sponse. 

We investigate the efficacy of ZT as a predictor of 
nonequilibrium device performance at the nanoscale by 
calculating the thermodynamic efficiency and power of 
an interacting quantum system using both nonequilib- 
rium many-body^ and Huckel theories. We discover that 
in both theories, variations of ZT and thermodynamic 
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FIG. 1. Near a 2n* order supernode in a device's transmis- 
sion spectrum, we find an order-dependent enhancement of 
the thermoelectric response which is limited only by the elec- 
tronic coherence length. Calculations were performed for a 
polyphenyl ether (PPE) SMJ with n repeated phenyl groups 
at room temperature (T=300K) with F=0.5eV. Notice that 
the enhancement is super-linear in n. Inset: ZT as a function 
jjb for n=l ... 5. 
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efficiency are in good qualitative agreement. However, 
large discrepancies between thermoelectric effects calcu- 
lated within many-body and Hiickel theory are found in 
the resonant tunneling regime, indicating the essential 
role of electron-electron interactions in nanoscale ther- 
moelectricity. For a thermoelectric quantum tunneling 
device, we find that the power output can be changed 
significantly by varying an external parameter, such as 
a gate voltage, and that this variation is not correlated 
with the variation of ZT . 

Neglecting inelastic processes, which are strongly sup- 
pressed at room temperature in SMJs, the current flowing 
into lead 1 of a two-terminal junction may be written as 
follows £ 



h ~h 



dE{E-yLiYT{E)[f 2 {E)-h{E)], (2) 



where v = (v=l) for the number (heat) current, f a (E) 
is the Fermi function for lead a with chemical potential 
fi a and inverse temperature j3 a , and T{E) is the trans- 
mission probability for an electron of energy E to tunnel 
across the junction. This transmission function may be 
expressed in terms of the junction's Green's functions 



as: 



T(E) = Tr {T 1 (E)G(E)T 2 (E)G' ! (E)} , 



(3) 



where T a (E) is the tunneling-width matrix for lead a 
and G(E) is the retarded Green's function of the SMJ. 

In organic molecules, such as those considered here, 
electron-phonon coupling is weak, allowing ZT to be ex- 
pressed as follows: 



ZT = ZT\ 



when 



.10 



ZT\ 



+ K ph / K el 
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and 



£M ( M) T)= J dE(E - T(E) 



(4) 



(5) 



(6) 



Here fo is the equilibrium Fermi function and 
K ph =KoT ph is the phonon's thermal conductance, 
where kq=(tt 2 /3)(k 2 i T/h) is the thermal conductance 
quantumii and T ph is the phonon transmission proba- 
bility. Since the Debye frequency in the metal lead is 
typically smaller than the lowest vibrational mode of a 
small organic molecule, the spectral overlap of phonon 
modes between the two is small, implying T ph< ^l and 
consequently that ZT^ZT\ el . 

Thermodynamically, a system's response is character- 
ized by the efficiency r\ with which heat can be converted 
into usable power V and the amount of power that can 
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FIG. 2. Schematic diagram of a thermoelectric device, where 
la is the heat current flowing into lead a, T a is the temper- 
ature and V is the power output. 

be generated. Applying the first law of thermodynamics 
to the device shown in Fig. ([2]) gives: 



(i) _ AO) 



= Ij. (Mi - M2), 



(7) 



where we mention that the power is equivalently phrased 
in terms of heat or electrical currents. The efficiency 
r\ is defined as the ratio of power output to input heat 
current: 



V 



r(i) 



r (i) 



n 



(i) 



(8) 



where we have assumed that Ti > T 2 . With these ex- 
pressions for the power and efficiency, we can completely 
quantify the performance of a quantum device, both near 
and far from equilibrium. 

As a first example, we calculate the non-linear ther- 
modynamic response of a meta-connected Au-benzcnc- 
Au SMJ using many-body 9 and Hiickel theory, shown 
in Fig. (|3aj) and Fig. (|3b|) . respectively. Although the 
transmission spectrum of this junction doesn't possess a 
supernode, it does possess a quadratic node within ir- 
electron theory^^ and will allow us to ascertain the im- 
portance of interactions on the thermoelectric response 
of a SMJ. 

In the top panel of each figure is a section of the trans- 
mission spectrum, showing the homo and LUMO reso- 
nances and the quadratic node directly in between at 
/U=/xo- Associated with this node is an enhancement 
in many linear-response metrics^ including ZT, which 
is shown in the second panel from the top. The bot- 
tom two portions of each figure show the calculated ef- 
ficiency 77 and power V when a junction with Xi=300K 
and T2=250K is further pushed out of equilibrium via 
the application of a bias voltage AV. In all simulations 
presented here, the lead-molecule coupling is taken to be 
symmetric such that V'^ lm =T5 na 5 ma , where n, m, and 
a are 7r-orbital labels and a is coupled to lead a. The 
efficiency is normalized with respect to the maximum al- 
lowed by the second law of thermodynamics, the Carnot 
efficiency 77c = AT/Ti, where AT=Ti-T 2 . 

The nonequilibrium thermodynamic response of a 1,3- 
benzene SMJ calculated using many-body theory is 
shown in Fig. (|3"a)) . The ZT and 77 spectra, shown in 
two middle panels of the same figure, exhibit peaks in 
the vicinity of both transmission nodes and resonances 
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FIG. 3. The transmission probability T(E), figure-of-merit ZT, Carnot-normalized efficiency r/frjc, and electrical power output 
V of a two terminal 1,3-benzene SMJ, with lead temperatures Ti=300K and T2=250K, calculated using (a) many-body and 
(b) Hiickel theory, highlighting the discrepancies near resonances and the similarities near the node in the two theories. As a 
function of [i, rj and ZT are in excellent qualitative agreement while V is only peaked near resonance, suggesting that ZT is 
incomplete as a device performance metric, (a) Many-body calculations give ■p poa i t =33/iW and »7 pe akA?C=ll-5% near resonance, 
(b) Hiickel calculations give 7 ? pea k=21/iW and r? pea k/^c=2.7% near resonance. The mid-gap region is discussed in Fig. Q. 
Note that the peak ZT— 0.75 is on par with currently available commercial thermoelectricsiA^ Calculations were performed 
using the model and parameterization of benzene discussed in detail in Ref.7 with F=0.63eV. 



whereas the power V , shown in the bottom panel, is only 
peaked near transmission resonances. Around either the 
HOMO or LUMO resonance, the peak power 7 ? pC ak=33^W 
and peak efficiency ^ poa k/ 7 ?c = ll-5% are only realized 
when the junction operates out of equilibrium at a bias 
voltage AF=3mV. With a chemical potential near the 
mid-gap node and AV=3.6mV 7? pe ak/??c=14.9%, larger 
than near resonance but with a much lower peak power 
^ P cak=.088nW. 

In the vicinity of a resonance, there are both quanti- 
tative and qualitative differences in the linear and non- 
linear thermodynamic response predicted by the two the- 
ories. By neglecting interactions, the Hiickel theory fails 
to accurately predict both the degeneracy and position 
of electronic resonances. It also incorrectly determines 
the peak values of ZT ', rj and V in the vicinity of a reso- 
nance. As can be seen near either (homo or lumo) res- 
onance in Fig. ©, the Hiickel theory predicts a Carnot- 
normalized peak efficiency of 2.7% which is nearly five 
times less than the 11.5% predicted by the many-body 
theory. The peak power near a resonance also varies 
considerably between the two theories, where the Hiickel 
calculations give 7 7 pC ak=21/iW while many-body theory 
predicts 7> poa k=33/jW. These results indicate that inter- 
actions are required to accurately predict the thermo- 
electric response of devices operating in the resonant- 
tunneling regime. It is interesting to note, however, that 
in both models the linear-response metric ZT qualita- 



tively captures the features of the non-linear metric rj. 

In this article, we are interested in thermoelectric en- 
hancement near nodes far away from any resonances. Al- 
though interactions are required in order to ensure the 
invariance of transport quantities under a global voltage 
shift (i.e. gauge- invariance), near the particle- hole sym- 
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FIG. 4. Calculations of ZT, rj and V in the vicinity of the 
transmission node at fj,—jio of a meta-benzene SMJ using 
many-body (red line and panel i) and Hiickel (black line and 
panel ii) theories, (a) and (b): ZT and rj are found to be 
identical and independent of theory, (c) V is strongly affected 
by interactions where, at peak efficiency (%,eak/»?c=14.91%), 
many-body and Hiickel calculations give 7 3 m ax=.088nW and 
7 ? max=1.87nW, respectively. The simulation parameters and 
colorscale are the same as in Fig. @. 



metric point the effect of interactions on the thermoelec- 
tric response should be small. In panels a-b of Fig. 
a comparison of ZT and r\ using both many-body and 
Hiickel theories is shown near fig for a 1,3-benzene SMJ. 
Near this point, ZT and 77 are independent of theory em- 
ployed. In contrast, the power, shown in panel c of the 
same figure, exhibits an order of magnitude difference be- 
tween the two theories. This observation can be under- 
stood by noticing that the calculated HOMO-LUMO gap 
is «10eV using many-body theory (panel c-i) whereas 
it is only «5.5eV when interactions are neglected in the 
Hiickel theory (panel c-ii). Since the power is peaked 
near transmission resonances, whose widths are fixed by 
the lead-molecule coupling T, the larger gap found us- 
ing many-body theory gives a correspondingly lower pre- 
dicted power. 

While the Hiickel theory is not able to accurately char- 
acterize the thermoelectric response of a junction in the 
resonant-tunneling regime, it is sufficient for predicting 
77 and ZT in the vicinity of the transmission node. Since 
we are interested in these quantities for mid-gap supern- 
odes, we shall use Hiickel theory to simulate the larger 
molecules presented below. 

The transmission node in a meta-benzene junction can 
be understood in terms of destructive interference of elec- 
tron waves traversing the ring at the Fermi energyJ^ Ac- 
cording to Luttinger's theorem , 13 ' 14 the Fermi volume is 
unaffected by the inclusion of electron-electron interac- 
tions. Consequently, in an aromatic ring such as benzene 
the Fermi wavevector fcp=7r/2d is conserved and is there- 
fore sufficient to characterize quantum interference both 
with and without interactions near since A(f>=kpAl, 
where A<f> is the relative phase between transport paths 
with length difference Al, and d is the inter-site distance. 

This is an important result, since the energy of res- 
onant levels will generally depend strongly on whether 
or not interactions are included. Since fcp is protected, 
however, the transmission node across a single phenyl 
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FIG. 5. A closeup of ZT and 77 near the quartic su- 
pernode of a 3,3'-biphenyl SMJ showing ZT pca k=1.84 and 
J7pcak/?7C=26.86% at a predicted power of 0.75pW. The junc- 
tion geometry is shown schematically in the inset of the upper 
panel. Simulations were performed using Hiickel theory with 
Ti=300K, T 2 =250K and T=0.5eV. 
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FIG. 6. Supernode enhancement of ZT, thermopower S and 
Lorenz number L for polyphenyl ether (PPE) SMJs with n 
repeated phenyl groups, shown schematically above the top 
panel. As a function of n, ZT pca t scales super-linearly ex- 
hibiting a peak value of 6.86 for n—6. The thermopower 
and Lorenz number are also enhanced with Spcaii—9571-iV/K 
and L pca k=55.33LwF at the same value of n. Simulations 
were performed using Hiickel theory at room temperature 
(T=300K) with F=0.5eV. Inter-phenyl electronic hopping was 
set an order of magnitude below the intra-phenyl value of 
2.64eV. 



group is not so much a coincidence of energy levels 
as a wave phenomenon, meaning that interference in 
molecules composed of multiple aromatic rings in series 
can be understood in terms of the interference within 
each subunit rather than the energy spectrum of the en- 
tire molecule. We find that such polycyclic molecules can 
exhibit higher-order supernodes, and that associated with 
a supernode is an order-dependent quantum enhance- 
ment of the junction's thermoelectric response. Addi- 
tional transport channels (e.g. cr-orbitals) or incoherent 
processes may lift the supernode. The effect on the ther- 
moelectric response is small provided the processes are 
weak, as discussed in Ref. (0). 

The 3,3'-biphcnyl junction, drawn schematically in the 
top panel of Fig. J5]), can viewed as two meta-connected 
benzene rings in series. This junction geometry is sim- 
ilar to that studied by Mayor et alJ^ In agreement 
with the prediction that a biphenyl junction should pos- 
sess a quartic supernode, the linear and non-linear re- 
sponse shown in Fig. ([5]) exhibits peak values of effi- 
ciency (77/770=26.86%) and ZT (1.84) that are over twice 
those of benzene. With ZTk-2, the biphenyl junction ex- 
hibits sufficient thermoelectric performance to be attrac- 
tive for many commerical solid-state heating and cooling 
applications^'^ As we shall see, this is only the first in 
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an entire class of supernode-possesing molecules which 
exhibit even larger values of r\ and ZT . 

In larger molecules composed of n meta-connected 
phenyl group in series, we expect that the transmission 
nodes should combine and give rise to a 2n th order su- 
pernode. Polyphenyl ether (PPE), shown schematically 
at the top of Fig. ^ , consists of n phenyl rings connected 
in series with ether linkages. Based on our previous dis- 
cussion, we predict that a PPE-based junction should 
exhibit a 2n th order supernode. The figure-of-merit ZT, 
thermopower S and Lorenz number L—n/GT for PPE 
junctions are shown in the top, middle and bottom pan- 
els of Fig. (J6j) , respectively, where the Lorenz number is 
normalized with respect to the Wiedemann-Franz (WF) 
value Lwf=(ti" 2 /3) (fce/e) 2 . 

The bottom panel of Fig. ^ shows an increasing 
peak Lorenz number L p0 ak with increasing n. In linear- 
response, L and S can be expressed in terms of Eq. ^ 
as: 



L \el ( eT )2 



£(0) 



£« 
£(0) 



(9) 



and S-— -^pjijoj, where e is the magnitude of the elec- 
tron's charge and T is the temperature. Using Eq. © 
and Eq. ([6]) with the transmission function of Eq. ([I) we 
find that: 



L 



WF 



[d* n+2 biTCSc(bir)] 

r2 



6=0 



[<9 b 2n 67rcsc(67r) 



(10) 



Setting n=6 in Eq. dTUJ) 

gives L max — 55.33-L\vF) corre- 



sponding exactly to the result of the full simulation 
shown in the bottom panel of Fig. ([6]). Similar agree- 
ment is found for the other values of n, confirming the 
presence of 2n th order supernodes in these junctions. 



We find that higher-order supernodes in the transmis- 
sion spectrum of a nanoscale junction give rise to an 
order-dependent quantum-enhancement of the linear and 
non-linear thermoelectric response. The full nonequilib- 
rium spectrum of thermodynamic efficiency qualitatively 
resembles the figure-of-merit ZT spectrum, suggesting 
that ZT encapsulates the salient physics related to effi- 
ciency even at the nanoscale. Efficiency, however, is only 
part of a device's performance. Another important quan- 
tity is the usable power produced by a device, whose vari- 
ations are poorly characterized by ZT at the nanoscale. 



Thermoelectric devices based on individual SMJs are 
ideally suited for local cooling in integrated nanoscale 
circuit architectures. Supernode-based devices have a 
low transmission probability and thus a large electri- 
cal impedance capable of withstanding voltage surges. 
Moreover, high-power macroscopic devices could be con- 
structed by growing layers of densely packed molecules. 
For example, a self-assembled monolayer with a sur- 
face density^ of 4x 10 15 molecules/cm 2 would give 
352kW/cm 2 at peak efficiency for a meta-benzene film. 
The efficiency of PPE-based devices increases with ring 
number and is only limited by the electronic coherence 
length, suggesting that highly efficient molecular-based 
thermoelectric devices may soon be realized. 
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